Quantum thermal equilibration from equipartition 
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The problem of mutual equilibration between two finite, identical quantum systems, A and B, 
prepared initially at different temperatures is elucidated. We show that the process of energy 
exchange between the two systems leads to accurate equipartition within energy shells in the Hilbert 
space of the total non-interacting, composite system, AiSB. This scenario occurs under the general 
condition of a weak interaction between the systems. We predict that the sole hypothesis of such 
equipartition is sufficient to obtain a relaxation of the peers, A and B, towards a common thermal- 
like state. This conjecture is fully corroborated by an exact diagonalization of several quantum 
models. 

PACS numbers: 05.30.-d,03.65.Aa 



I. INTRODUCTION 

The time evolution of an isolated quantum system af- 
ter applying a sudden change for one of its parameters, 
i.e., - a quench - has recently gained considerably at- 
tention, both in the theoretical and experimental physics 
communities [l[. State of art numerical simulations Q- 
|8|, motivated by recent advances in manipulations with 
ultracold atoms Q, have not only allowed to validate a 
number of theoretical predictions pi]| - [l2| . but also pro- 
duced several conceptually new research directions. One 
of these tracks refers to the exploration of the quench 
machinery as an effective tool to drag the system of in- 
terest into a new state. The latter can effectively mimic 
the state of thermal equilibrium - without the need of 
coupling the system to a heat bath [IJ]. 'Mimic' means 
here that the expectation values of relevant observables 
are close to those following from the thermal Gibbs state, 
QT(xexp{-l3H),p = l/kBT. 

The equilibration between two identical, initially non- 
interacting systems, A and B, can be considered as a 
quench applied to the composite system, 



H^ = HA(E)lB + 'i-A®HB + X{t)H" 



starting out from the noninteracting limit, A = 0, to the 
regime of interaction, A = Aint- It has been shown with 
prior work [l4| that for the initial product state, prepared 
at different temperatures, Ta and Tb, g{0) = Q^ ® Q^ , 
the step-like quench \{t) — Aint^(i) evolved the com- 
posite system into a new state, Q{t), such that, for the 
times t > ieq, the reduced density matrices, g^{t) and 
g^{t), become quasistationary [l5| and mimic perfectly 
a thermal equilibrium with a common temperature Toq. 
Although this scenario seemingly is universal, in a sense 
that it works equally well for very different physical sys- 
tems, the physical mechanism at work remained elusive. 
With this study we address this open problem. We 
show that mutual thermal relaxation of two finite quan- 
tum systems follows from a generic hypothesis about the 
asymptotic state of the composite system after applica- 
tion of a weak interaction quench: namely, equipartition 



inside energy shells E = e + e^ of the identical spectra 
of the composite system Ha ® Hb constitutes a sufficient 
condition for the emergence of the mutual thermal equi- 
libration between the system's halves. We corroborate 
this conjecture by using four different types of models, 
including synthesized Hamiltonians with different distri- 
butions of energy levels and a system of two interacting 
spin clusters. 



II. SETUP 

The model ([T]) consists two identical quantum systems, 
A and B, with identical finite spectra, {cj}, j = 0, ..., A^ — 
1, of width Ae = e^v-i — eg, and a set of eigenstates {|j)}. 
The corresponding energy level distribution is encoded by 
the density of states. 



N-\ 



n(e) = ^ 5(e - ej). 



(2) 



j=o 



The initial states of the systems are given by Gibbs den- 
(1) sity matrices, Q^ifi) = Qta ^-^d Q^{Q) — QTb^ ^-t the tem- 
peratures Ta and Tb- The initial state of the total com- 
posite system in the product basis, {|^^„(j ■•)) = \i)®\k)}, 



E. 



m{j,k) 



efc, is represented by a diagonal den- 



A „B _ X. ..X. ..^0 
jk 

Z7 Z^ , where Za b are the 



PTk\]'k' = PT,j'Pk,k' = ^j,3'^k,k'P"jk 



sity matrix, p*°*(0) 

P% = exp{-ejl3A - ekl3B)^A -b 
partition functions, Za.b = /_ dec'xp{—e/3A,B)'n{e). 
Henceforth, we set eg = and use Ae as the energy unit 
if not specified otherwise. 



III. EQUILIBRATION INDUCED BY 
EQUIPARTITION 

We define the energy shells of the composite system 
in the product basis {!*„)} by using the condition 
\Em — E\ < S [ial- The constant 5 is chosen small with 
respect to the spectral width, S <C Ae, but still larger 



than the mean level spacing of the composite system, 
S > 2Ae/(7V - 1), with M = N x N energy levels. 
The last condition implies that the every energy shell 
contains many eigenstates. The switch-on of an inter- 
action Hamiltonian, Aint-ff'"*, which is non-diagonal in 
the product basis, generates a set of new eigenstates, 
|*m): ^^l*m) = -E^|*m)- If we sort both sets of 
eigenstates, {1^^)} ^^^ il'^m)}, with respect to their 
energies, Em and E^, we obtain a bell-shaped overla 
function /m'(m) = |('J'^'l^m'+'ra)P' centered at m 
with a width that grows with the strength of perturba- 
tion Aint I13i im ■ Throughout this study we assume the 
weak coupling limit. 
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C) « Ae , 



(3) 



obeying, in addition, the condition 

Ai„t||i/'"V'°'(0)|| > stot - 2Ae/(AA- 1), (4) 

where ||...|| is the operator norm in the Hilbert space of 
the composite system and Stot is the mean level-spacing. 
The last condition means that the interaction should 
not be too weak, otherwise the non-thermal scenario of 
arithmetic-mean equilibration [14| would take place. 

In common setups of quench studies the isolated sys- 
tem is initially prepared in a m-th eigenstate (typically in 
its ground state, |^[]) [i-Ii,[il) of the Hamiltonian H^=°. 
A weak quench then results in a local smearing of the ini- 
tial wave function over the narrow set of new eigenstates, 
given by the function fm'{m), so that 'microcanonical 
thermalization' can be expected 0, [H, [ill, [H H^] • Mi- 
crocanonical thermalization implies that a closed quan- 
tum system is transformed into a new state, which sa tisfy 
Boltzmann's postulate of equal a priori probability |16l |. 
here applied to the quantum states within an energy 
shell [21|. Evidence that this may indeed be expected 
under quite general conditions [^, [2J, [2J] is nowadays 
discussed under the label "quantum typicality" [2^, |29| ■ 

We start by extending the concept of thermalization to 
the case of the bipartite system initially prepared in the 
product state g{0). At time i = we turn on the quench 
by setting A ^ 0. Then, after some elapsed characteris- 
tic time iroi, we switch-off the perturbation and investi- 
gate the state of the system with respect to the product 
basis |^°j). By representing the system Hilbert space, 
sheared by energy shells of different energies E, as having 
an onion-like structure, we conjecture that a proper per- 
turbation will initiate the population exchange between 
the eigenstates within each shell, - independently of the 
remaining part of the system Hilbert space 



271. This 



exchange will lead finally to the equipartition of the level 
populations within each energy shell. 

In order to cast our hypothesis into a formal math- 
ematical language we first introduce a two-dimensional 
probability density function (pdf): 



j,k=0 



6,)<5(6^-6,), (5) 




After equilibration 
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FIG. 1: (color online). Thermal equilibration between two fi- 
nite, identical quantum systems following the equipartition 
scenario. Systems, initially prepared at different temper- 
atures (left), eventually arrive after relaxation at a quasi- 
equilibrium state (right), characterized by distributions uni- 
form along the E-axis (equipartition). The color coding (de- 
creasing in weight from bright to dark) depicts the behavior of 
the pdf P[e^, e^](f), see Eq. ^. Note that the equipotential 
lines on the left panel are not horizontal but slightly inclined. 
The inclination is small due to a large difference between the 
peer's initial temperatures, Ta 3> Tb- 



where the populations Pjk(t) are governed by the di- 
agonal elements of the total density matrix P^Ikj'k'W- 
The initial pdf is given by P[ef,ef]{0) = exp(-/JAef - 
(3B£f)ZA^ZB^, see Fig. [T] It is useful to introduce the 
auxiliary variables, E = e^ + e^ and T, = e^ ~ e^ , which 
form a new coordinate axes. The first variable, E, de- 
fines the above mentioned energy shell, while the second 
one, E, can be used to label the states within the shell. In 
this representation the initial condition assumes the form 
pO^ = P0(£;™,E„) = cxp{-Eml3+ - T.mP-)Z^^Zg\ 
with the two inverse temperatures /?+ = {j3a + /3b)/2 
and /3^ = (/3b — /3/i)/2. The density of states in new 
variables, fi{E,iy), does generally not factorize. 

According to the proposed equipartition scenario, after 
equilibration the diagonal elements of the total system 
density matrix, P'"i{E,T,) = P'=<i(i;), derive from the 
equipartition of the probability over the corresponding 
energy shells, reading 



p'=i(£;) 



-^-^ rj;'^ljn{E,i:)e-p-^di: 



ZaZ 



A^B 



r_}^%^n{E,i:)d^ 



(6) 



where the integration limits are ri{E) ^ E ior < E < 
Ae, and t]{E) = 2Ae - E ioi Ae < E < 2Ae, see Fig. 1. 
Note that the expression (j6]) conserves energy within a 
specific shell. Then the energy level populations of a 
single peer can be evaluated as: 



Pi 



P'^'^iej + e)n{e)de. 



(7) 



Note that the equilibrium distribution explicitly depends 
on the density of states, n{e), of the system Hamiltonian. 
Below, by using three different classes of system Hamil- 
tonians we demonstrate that (i) the smearing along the 
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FIG. 2: (color online). Equilibration between two identical 
finite quantum systems with A'^ = 151 uniformly distributed 
energy levels. In panel (a) we depict the diagonal elements of 
the total system density matrix before the interaction quench 
and in panel (b) the result after equilibration occurred. Notice 
full agreement with the equipartition scenario, see Fig. 1. The 
initial temperatures are keTA = 13 a ~ 2Ae and kuTs ~ 
Pg^ = O.lAe. Panel (c) depicts the populations before and 
after equilibration. The initial populations, Pi (0) and Pi (0), 
are denoted by the thick solid (blue for system A and red for 
system B) lines. The resulting equilibrium populations (thick 
dashed line) agree (within line thickness) with the analytical 
prediction, Eq. Q, and are very close to the canonical thermal 
populations obtained from the energy conservation condition, 
Eq. ((H]) (thin (green) line). Energy is measured in units of the 
spectral width Ae. 



E-axis is sufficient in producing thermal equilibration be- 
tween the peers and in fact (ii) such smearing indeed is 
achieved in those systems after an interaction quench. 

In order to validate our predictions we performed cal- 
culations for three different classes of synthesized Haniil- 
tonians, with uniform, semicircular and a Gaussian den- 
sity of states, Eq. ^. Finally, we investigated the ther- 
mal equilibration between two finite spin clusters. 



IV. THERMAL EQUILIBRATION BETWEEN 

PEERS WITH UNIFORM DISTRIBUTIONS OF 

ENERGY LEVELS 



within the weak-coupling limit, so that the interaction 
quench does not cause appreciable heating of the com- 
posite system, but still is strong enough as to guarantee 
the thermal-like equilibration scenario [IJ]. Here we use 
the dimensionless coupling constant Aint = 0.015(s/ft.), 
where s = Ae/{N — 1) is the mean level spacing, and 

The initial and the equilibrium population pdf 's, as ob- 
tained by the exact diagonalization of the composite sys- 
tem with Af = N X N ~ 22801 states, are presented with 
Fig. m The equilibrium pdf P{E, S) shown in Fig. ^h) 
assumes a stripe-like structure, being uniform along the 
S-axis, in full agreement with the equipartition scenario, 
see Fig.[T] We also checked that the emerging equilibrium 
populations P'^'^{E) follow closely the predicted result in 

Eq. dg. 

The equilibrium values oip°'^ for a single peer obtained 
by using Eq. ([7]) are shown in Fig. [2][c) by the dashed 
thick line. The analytical prediction in Eq. ([7]) are in- 
distinguishable from the numerical data points obtained 
from the direct diagonalization of the composite Hamil- 
tonian in Eq. ([T]) . Except for some small deviation in the 
high-energy tail, both distributions fit almost perfectly 
the thermal distribution with the equilibrium tempera- 
ture Ti^ri extracted from the condition of energy conser- 
vation Q, 






^cq 



g — efc/fc^T^ g- 

ffc — ?77? hefe — 



k/keTB 



2Za 



2Zf 



,(8) 



see the thin solid (green) line in Fig. [JJc). 



V. SEMICIRCULAR DISTRIBUTION OF 
ENERGY LEVELS 

In the present example we use spectra that are typ- 
ical for the class of Hamiltonians modeled by a ran- 
dom matrix drawn from GOE [30[ . In the limit of large 
number of levels, the density of states of a single peer 
can be approximated by the continuous semicircular dis- 
tribution, n(e) = (4/7r) a/1/4 - (e - 1)^/4, see the in- 
set in Fig. ^a.). Thus, for the total system we have 
n{E,T,) = l/ny/{E'^ - £2)2 - A{E^ - Y.^)[E - 1). The 
results of the exact diagonalization perfectly match the 
prediction in Eq. ([T]). As for the first example, the thin 
solid (green) line indicates the distribution with the equi- 
librium temperature given by Eq. ([5]). 



We have synthesized a Hamiltonian with A^ == 151 lev- 
els, distributed them randomly and uniformly in the in- 
terval [0, Ae]. The interaction Hamiltonian is composed 
as the product of two identical matrices, iJ'"* = Ya (g) Yg, 
where the matrix Y = Ya ^ Yb has been drawn from 
a Gaussian Orthogonal Ensemble (GOE). Namely, Y = 
{R + R^)/2, where the matrix R in the product basis 
is given by its real elements obeying standard normal 
distribution [2^. The interaction between the peers is 



VI. MODEL WITH GAUSSIAN DISTRIBUTION 
OF ENERGY LEVELS 

This next class of Hamiltonians refers to quantum sys- 
tems possessing a finite number of interacting particles 
or spins, as realized with fermionic 31[ and bosonic Hub- 
bard models [S^]. In the limit N — > cxd the correspond- 
ing density of states can be approximated by the contin- 




FIG. 3: (color online). The same as in Fig. [2] but here for the 
equilibration scenario between two identical peers with a semi- 
circle (a) and Gaussian (b) density of states. The variance of 
the Gaussian distribution is cr = 0.1833Ae. The remaining 
parameters are the same as in Fig. [Jfc). Insets: exact semi- 
circle and Gaussian distributions, dashed (black) lines, and 
the density of states of the finite synthesized Hamiltonian 
with A'^ = 181 energy levels (histograms). 



uous Gaussian function ?i{e) oc cxp [— (e — l/2)^/(2o'-^)] , 
wherein both the width a and energies e are in units of 
the total width Ae. 

In distinct contrast to the semicircle distribu- 
tion, the Gaussian density of states remains fac- 
torized after the frame transformation, n(E, S) ex 
exp [^{E - 1/2)2/(20-2)] exp [-E^/{2ct^)] . Therefore, 
Eq. ([5]) reduces (up to irrelevant normalization constant) 
to the form: 



P"'i(£:) ex e" 



iJ-viE)^ a^ 



r,{E) 



J-t; 



-SV(2a2)dS 



(9) 



In the limit of a very broad Gaussian distribution, a ^ 
Ae, the above expression approaches the foregoing result 
of a uniform distribution, see Fig. [2l In the opposite 
limit of a very narrow distribution; i.e., a <C ^Ae//3~, 
the integrals in the numerator and denominator of Eq. ^ 
yield approximately the same values, thus rendering the 
Boltzmann-like distribution. 



P'^'^iE) ex e 



-P+E 



(10) 



for the composite system. This limit corresponds to a 
"strong thcrmalization" numerically observed with two 
coupled Bose- Hubbard models [SJ]. Accordingly, both 
peers also relax to the thermal states of the same tem- 
perature, Teq = (fcB/3+)~^, see Fig.|3|b). It is noteworthy 
that the strong thcrmalization was absent in the previ- 
ously considered cases. 

To conclude this section, we discuss the important is- 
sue of off-diagonal elements of the reduced density ma- 
trices of the peers after they reached the state of a joint 




FIG. 4: (color online). The reduced density matrices of peers 
with the uniform density of states after equilibration is com- 
pleted. Crosses (red) show the elements of the density nratrix 
of system A (initially "hot" ) , while open circles (blue) show 
those for system B (initially "cold"). Both peers were initially 
in canonical thermal states, with the corresponding density 
matrices in diagonal form. The parameters are the same as 
those in Fig. [D 



thermal equilibrium. With Figs. [T][3l we addressed the 
diagonal elements of the reduced density matrices only 
and showed that they fit the thermal distributions with 
the equilibrium temperature given by Eq. (|5]). Remark- 
ably, the off-diagonal elements, although they appear 
during the equilibration process, remain extremely small 
after equilibration is completed. Therefore, the thermal- 
ized density matrices of the peers preserve their diagonal 
forms and remain near canonical, see Fig. |4l 



VII. THERMAL EQUILIBRATION OF 
INTERACTING SPIN CLUSTERS 

Synthesized Hamiltonians, although very useful for nu- 
merical studies 



27[, have a serious drawback. Namely, 



they do not feature some nontrivial statistical properties 
which may be present in spectra of actual quantum sys- 
tems. Therefore the equipartition scenario needs to be 
tested with a realistic physical Hamiltonian. 

As a last peer model we use a finite cluster of Ns = 
8 interacting 1/2-spins. Two clusters are placed into a 
constant magnetic field, pointing along z-direction, and 
brought into a local contact, see Fig. [5ja). Each cluster 
has A^ = 2^ = 256 states, so that the overall dimension of 
the Hilbcrt space of the composite system is TV = 2^^^ ~ 
216 = 65536. 

For two identical clusters, we employ here the spin 
model that is also referred as to XXZ model with the 
following Hamiltonian, Ha = Hb = H: 

H = vJ2stS]-jY,{S^S^ + SfSf) + MY,St, (11) 
where Sf,Sf,S^ are spin-1/2 operators on site i, V (J) 



(a) 



Spin cluster A 




FIG. 5: (color online), (a) Two locally interacting spin clus- 
ters, (b) The diagonal elements of the density matrices cor- 
responding to the different magnetization subspaces after the 
thermalization process is completed. The Hilbert space of a 
single cluster splits into 2Ns + 1 = 9 invariant subspaces, 
"Hs^ , X — A,B, with a spin Sx taking integer values from 
—Ns/2 = —4 to Ns/2 = 4. The interaction between clusters 
leaves the total magnetization of the composite system, Sz, 
invariant. Therefore, for the chosen initial conditions with 
subspaces S% s = only populated, the consecutive equili- 
bration process is restricted to the Sz = subspace of the 
Hilbert space of the composite system. As a result of the in- 
teraction, all possible products of local subspaces, Tis^ ^Tis' , 
with opposite magnetization, S\ = —S%, become populated. 
Note that the single-state subspaces with S^'^ = ±4 are not 
shown. Initial temperatures of the clusters are the same as in 
Fig.H 
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FIG. 6: (color online). Populations of different S% g sub- 
spaces of the Hilbert space of a single cluster after the equi- 
libration process is completed. Top panels: Equilibrium pop- 
ulations of the energy levels in subspaces S%, being marked 
by the (blue) crosses, and S% — —S\, as marked by (red) 
circles, are shown for the cluster A and for the cluster B, re- 
spectively. The solid line for each subspace depicts the ther- 
mal energy level populations at the equilibrium temperature 
Teq, multiplied by the total population of the corresponding 
subspace. The latter subspace is marked by the filled grey 
symbol in each panel, as shown in the bottom part. The 
equilibrium populations and the thermal distributions agree 
(within line thickness) with the analytical prediction. Bot- 
tom panel: Individual population values of the corresponding 
subspaces S% — —S% after equilibration. The parameters are 
the same as in Fig. (5] 



are the exchange constants in z {x, y) directions, M is the 
external magnetic field, and (...) indicates here all pairs 
of next-neighbor spins connected according to bonds of 
a single spin cluster displayed in Fig. [SJa). The coupling 
term between the clusters, iJ™*, assumes similar to the 
Hamiltonian of the spin cluster form: 



ij"t = vY. sts^ -jJ2 (sfsj + sfsy). (12) 



E 



to);. 



Here the sum runs over the two bonds {ij)x that bind 
the two spin clusters together upon the action of quench. 

In distinct contrast to the synthesized model discussed 
before, both single clusters and the entire composite sys- 
tem possess integrals of motion additional to the total 
energy. That are the total magnetization along direction 
of the applied magnetic field, S'^ ^ = ^ja b ^h ^°'" ^^^ 
clusters, and z-component of the total spin, S"^ = ^ 5*1, 
for the composite system [2^. As a consequence, the 
Hamiltonian of a single X-cluster, X = Aoy B, factorizes 
over the product space T-Ls" into 2A^s -I- 1 independent 
blocks. So does the Hamiltonian of the composite system 
over the product space ^Hs^, yielding ANs + 1 blocks. 

Conservation of the total magnetization allows to 
study the process of mutual quantum equilibration in 



a more complex situation. For both clusters we choose 
initial states with only invariant subspaces S 



A,B 





thermally populated. By resorting to the equipartition 
hypothesis, we predict that a weak interaction quench 
that preserves the magnetization of the composite sys- 
tem, 5*^ = S% + Sg ~ 0, but violates the separate conser- 
vation of the magnetization of individual cluster, S\ g, 
would not only lead to the equilibration between the sub- 
spaces S\ B — 0, but shall also initiate a population and 
consecutive thermalization within subspaces S\ ^ 7^ 0. 

Our analytical calculations based on generalized 
form of Eqs. ([S][7]) for the factorized space ^{5^=0 = 
J2s' '^S' ® 'H-S' perfectly agree with exact diagonal- 
ization of the model Hamiltonian in the subspace of zero 
total magnetization, S^' ~ 0, spanned by 12870 states, 
see Figs. [5l [51 The model parameters are J = 0.2Ae, 
V = O.lAe, M = O.OSAe, Ai„t = 0.095(s//i) = 1. 

The equilibrium temperature T^q was calculated by us- 
ing Eq. ([5]), which was applied to the initially populated 
subspace, 'Hs'=o (SiHs'^o, only. It is noteworthy that 
the 'equilibrium' distributions for different subspaces per- 
fectly match the thermal distributions with the same 
equilibrium temperature, Toq, see in Fig. [H] (top panels). 



VIII. SUMMARY AND OUTLOOK 

In conclusion, using different classes of Hamiltonians, 
we have unraveled the mechanism responsible for the 
thermal equilibration of two identical quantum peers pre- 
pared initially in canonical states at different tempera- 
tures. This mechanism, i.e., the equipartition within en- 
ergy shells in the Hilbert space of the composite system, 
may appear whenever the interaction is small enough to 
satisfy the weak-coupling condition, given by Eqs. ([3j 
|4]). However, the equipartition scenario is not univer- 
sal: Quantum systems that exhibit Anderson localiza- 
tion are expected to invalidate the equipartition scenario 
when coupled by a weak local interaction, and the final 
equilibrium states of the correspondi ng p eers can differ 
substantially from being thermal-like [3^, |3Q| . 

One should keep in mind that the time evolution of any 
isolated quantum system with a finite number of levels 
has a finite recurrence time, which depends on the sys- 
tem spectrum and the system initial conditions. Thus 



the equilibration of the peers to a thermal 'equilibrium' 
after some interaction time t does not contradict the 
disappearance of the equilibration at some larger times, 
irec > t, due to revivals. The revival time scales can be 
very short when the interacting systems are small j37| . 

The equilibration process is governed by the Hamilto- 
nians, Ha, Hb, and _ff'"*, and its output is in one-to-one 
correspondence with the initial states of the peers. It 
means that initial states different from thermal Gibbs 
states, generally would lead to a final quasi-equilibrium 
which may not be thermal-like anymore. This compli- 
cation, however, could be weakened by the increase of 
the number of peers: interaction between M ^ 2 sys- 
tems would effectively mimic an environment for a single 
peer, thus leading to the mutual equilibration of all peers 
to nearly identical thermal states regardless the shape of 
their initial eigenstate distributions j39l |40[ . 

We acknowledge the support by the German Excel- 
lence Initiative "Nanosystems Initiative Munich (NIM)". 
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